Genome-wide association study (GWAS) reveals polygenic architecture for limber pine quantitative disease resistance to white pine blister rust.

Development of durable resistance effective against a broad range of pathotypes is crucial for restoration of pathogen-damaged ecosystems. This study dissected the complex genetic architecture for limber pine quantitative disease resistance (QDR) to Cronartium ribicola using a genome-wide association study (GWAS). Eighteen-month-old seedlings were inoculated for resistance screening under controlled conditions. Disease development was quantitatively assessed for QDR-related traits over four years post inoculation. To reveal genomic architecture contributing to QDR-related traits, a set of genes related to disease resistance with genome-wide distribution was selected for targeted sequencing for genotyping of single-nucleotide polymorphisms (SNPs). GWAS revealed a set of SNPs significantly associated with quantitative traits for limber pine QDR to WPBR, including number of needle spots and stem cankers, as well as survival four years post-inoculation. The peaks of marker-trait associations displayed a polygenic pattern with genomic regions as potential resistant quantitative trait loci (QTLs), distributed over ten of the 12 linkage groups (LGs) of Pinus. None of them were linked to the Cr4-controlled major gene resistance (MGR) previously mapped on LG08. Both normal canker and bole infection were mapped on LG05, and the associated SNPs explained their phenotypic variance up to 52%, tagging a major resistant QTL. Candidate genes containing phenotypically associated SNPs encoded putative nucleotide-binding site leucine-rich repeat (NBS-LRR) proteins, LRR-receptor-like kinase (LRR-RLK), cytochrome P450 superfamily protein, heat shock cognate protein 70, glutamate receptor, RNA-binding family protein, and unknown protein. The confirmation of resistant QTLs broadens the genetic pool of limber pine resistance germplasm for resistance breeding.

Transcriptomic Reprogramming and Genetic Variations Contribute to Western Hemlock Defense and Resistance Against Annosus Root and Butt Rot Disease.

Western hemlock (Tsuga heterophylla) is highly susceptible to Annosus root and butt rot disease, caused by Heterobasidion occidentale across its native range in western North America. Understanding molecular mechanisms of tree defense and dissecting genetic components underlying disease resistance will facilitate forest breeding and disease control management. The aim of this study was to profile host transcriptome reprogramming in response to pathogen infection using RNA-seq analysis. Inoculated seedlings were clearly grouped into three types: quantitative resistant (QR), susceptible (Sus), and un-infected (Uif), based on profiles of H. occidentale genes expressed in host tissues. Following de novo assembly of a western hemlock reference transcriptome with more than 33,000 expressed genes, the defensive transcriptome reprogramming was characterized and a set of differentially expressed genes (DEGs) were identified with gene ontology (GO) annotation. The QR seedlings showed controlled and coordinated molecular defenses against biotic stressors with enhanced biosynthesis of terpenoids, cinnamic acids, and other secondary metabolites. The Sus seedlings showed defense responses to abiotic stimuli with a few biological processes enhanced (such as DNA replication and cell wall organization), while others were suppressed (such as killing of cells of other organism). Furthermore, non-synonymous single nucleotide polymorphisms (ns-SNPs) of the defense- and resistance-related genes were characterized with high genetic variability. Both phylogenetic analysis and principal coordinate analysis (PCoA) revealed distinct evolutionary distances among the samples. The QR and Sus seedlings were well separated and grouped into different phylogenetic clades. This study provides initial insight into molecular defense and genetic components of western hemlock resistance against the Annosus root and butt rot disease. Identification of a large number of genes and their DNA variations with annotated functions in plant resistance and defense promotes the development of genomics-based breeding strategies for improved western hemlock resistance to H. occidentale.

Comparative Association Mapping Reveals Conservation of Major Gene Resistance to White Pine Blister Rust in Southwestern White Pine (Pinus strobiformis) and Limber Pine (P. flexilis).

All native North American white pines are highly susceptible to white pine blister rust (WPBR) caused by Cronartium ribicola. Understanding genomic diversity and molecular mechanisms underlying genetic resistance to WPBR remains one of the great challenges in improvement of white pines. To compare major gene resistance (MGR) present in two species, southwestern white pine (Pinus strobiformis) Cr3 and limber pine (P. flexilis) Cr4, we performed association analyses of Cr3-controlled resistant traits using single nucleotide polymorphism (SNP) assays designed with Cr4-linked polymorphic genes. We found that ∼70% of P. flexilis SNPs were transferable to P. strobiformis. Furthermore, several Cr4-linked SNPs were significantly associated with the Cr3-controlled traits in P. strobiformis families. The most significantly associated SNP (M326511_1126R) almost colocalized with Cr4 on the Pinus consensus linkage group 8, suggesting that Cr3 and Cr4 might be the same R locus, or have localizations very close to each other in the syntenic region of the P. strobiformis and P. flexilis genomes. M326511_1126R was identified as a nonsynonymous SNP, causing amino acid change (Val376Ile) in a putative pectin acetylesterase, with coding sequences identical between the two species. Moreover, top Cr3-associated SNPs were further developed as TaqMan genotyping assays, suggesting their usefulness as marker-assisted selection (MAS) tools to distinguish genotypes between quantitative resistance and MGR. This work demonstrates the successful transferability of SNP markers between two closely related white pine species in the hybrid zone, and the possibility for deployment of MAS tools to facilitate long-term WPBR management in P. strobiformis breeding and conservation.

Fine dissection of limber pine resistance to Cronartium ribicola using targeted sequencing of the NLR family.

BACKGROUND: Proteins with nucleotide binding site (NBS) and leucine-rich repeat (LRR) domains (NLR) make up one of most important resistance (R) families for plants to resist attacks from various pathogens and pests. The available transcriptomes of limber pine (Pinus flexilis) allow us to characterize NLR genes and related resistance gene analogs (RGAs) in host resistance against Cronartium ribicola, the causal fungal pathogen of white pine blister rust (WPBR) on five-needle pines throughout the world. We previously mapped a limber pine major gene locus (Cr4) that confers complete resistance to C. ribicola on the Pinus consensus linkage group 8 (LG-8). However, genetic distribution of NLR genes as well as their divergence between resistant and susceptible alleles are still unknown. RESULTS: To identify NLR genes at the Cr4 locus, the present study re-sequenced a total of 480 RGAs using targeted sequencing in a Cr4-segregated seed family. Following a call of single nucleotide polymorphisms (SNPs) and genetic mapping, a total of 541 SNPs from 155 genes were mapped across 12 LGs. Three putative NLR genes were newly mapped in the Cr4 region, including one that co-segregated with Cr4. The tight linkage of NLRs with Cr4-controlled phenotypes was further confirmed by bulked segregation analysis (BSA) using extreme-phenotype genome-wide association study (XP-GWAS) for significance test. Local tandem duplication in the Cr4 region was further supported by syntenic analysis using the sugar pine genome sequence. Significant gene divergences have been observed in the NLR family, revealing that diversifying selection pressures are relatively higher in local duplicated genes. Most genes showed similar expression patterns at low levels, but some were affected by genetic background related to disease resistance. Evidence from fine genetic dissection, evolutionary analysis, and expression profiling suggests that two NLR genes are the most promising candidates for Cr4 against WPBR. CONCLUSION: This study provides fundamental insights into genetic architecture of the Cr4 locus as well as a set of NLR variants for marker-assisted selection in limber pine breeding. Novel NLR genes were identified at the Cr4 locus and the Cr4 candidates will aid deployment of this R gene in combination with other major/minor genes in the limber pine breeding program.

Comparative Transcriptomics and RNA-Seq-Based Bulked Segregant Analysis Reveals Genomic Basis Underlying Cronartium ribicola vcr2 Virulence.

Breeding programs of five-needle pines have documented both major gene resistance (MGR) and quantitative disease resistance (QDR) to Cronartium ribicola (Cri), a non-native, invasive fungal pathogen causing white pine blister rust (WPBR). WPBR is one of the most deadly forest diseases in North America. However, Cri virulent pathotypes have evolved and can successfully infect and kill trees carrying resistance (R) genes, including vcr2 that overcomes MGR conferred by the western white pine (WWP, Pinus monticola) R gene (Cr2). In the absence of a reference genome, the present study generated a vcr2 reference transcriptome, consisting of about 20,000 transcripts with 1,014 being predicted to encode secreted proteins (SPs). Comparative profiling of transcriptomes and secretomes revealed vcr2 was significantly enriched for several gene ontology (GO) terms relating to oxidation-reduction processes and detoxification, suggesting that multiple molecular mechanisms contribute to pathogenicity of the vcr2 pathotype for its overcoming Cr2. RNA-seq-based bulked segregant analysis (BSR-Seq) revealed genome-wide DNA variations, including about 65,617 single nucleotide polymorphism (SNP) loci in 7,749 polymorphic genes shared by vcr2 and avirulent (Avcr2) pathotypes. An examination of the distribution of minor allele frequency (MAF) uncovered a high level of genomic divergence between vcr2 and Avcr2 pathotypes. By integration of extreme-phenotypic genome-wide association (XP-GWAS) analysis and allele frequency directional difference (AFDD) mapping, we identified a set of vcr2-associated SNPs within functional genes, involved in fungal virulence and other molecular functions. These included six SPs that were top candidate effectors with putative activities of reticuline oxidase, proteins with common in several fungal extracellular membrane (CFEM) domain or ferritin-like domain, polysaccharide lyase, rds1p-like stress responsive protein, and two Cri-specific proteins without annotation. Candidate effectors and vcr2-associated genes provide valuable resources for further deciphering molecular mechanisms of virulence and pathogenicity by functional analysis and the subsequent development of diagnostic tools for monitoring the virulence landscape in the WPBR pathosystems.

In-vitro anti-fungal assay and association analysis reveal a role for the Pinus monticola PR10 gene (PmPR10-3.1) in quantitative disease resistance to white pine blister rust.

Pathogenesis-related (PR) proteins play important roles in plant defense response. However, functional investigation of PR10 genes is still limited and their physiological roles have not been conclusively characterized in biological processes of conifer trees. Here, we identified multiple novel members in the western white pine (Pinus monticola) PmPR10 family by bioinformatic mining available transcriptomic data. Phylogenetic analysis of protein sequences revealed four PR10 and two PR10-like clusters with a high synteny across different species of five-needle pines. Of 10 PmPR10 genes, PmPR10-3.1 was selected and expressed in Escherichia coli. The purified recombinant protein exhibited inhibitory effects on spore hyphal growth of fungal pathogens Cronartium ribicola, Phoma exigua, and Phoma argillacea by in-vitro anti-fungal analysis. Genetic variation analysis detected a total of 21 single nucleotide polymorphisms (SNPs) within PmPR10-3.1 in a collection of P. monticola seed families. A nonsynonymous SNP (t178g) showed significant association with relative levels of quantitative disease resistance (QDR), explaining about 8.7% of phenotypic variation as the peak value across all SNPs. Our results provide valuable insight into the genetic architecture underlying P. monticola QDR and imply that PmPR10-3.1 may function as an important component in conifer basal immunity for non-specific resistance to a wide spectrum of pathogens.

Association Mapping and Development of Marker-Assisted Selection Tools for the Resistance to White Pine Blister Rust in the Alberta Limber Pine Populations.

Since its introduction to North America in the early 1900s, white pine blister rust (WPBR) caused by the fungal pathogen Cronartium ribicola has resulted in substantial economic losses and ecological damage to native North American five-needle pine species. The high susceptibility and mortality of these species, including limber pine (Pinus flexilis), creates an urgent need for the development and deployment of resistant germplasm to support recovery of impacted populations. Extensive screening for genetic resistance to WPBR has been underway for decades in some species but has only started recently in limber pine using seed families collected from wild parental trees in the USA and Canada. This study was conducted to characterize Alberta limber pine seed families for WPBR resistance and to develop reliable molecular tools for marker-assisted selection (MAS). Open-pollinated seed families were evaluated for host reaction following controlled infection using C. ribicola basidiospores. Phenotypic segregation for presence/absence of stem symptoms was observed in four seed families. The segregation ratios of these families were consistent with expression of major gene resistance (MGR) controlled by a dominant R locus. Based on linkage disequilibrium (LD)-based association mapping used to detect single nucleotide polymorphism (SNP) markers associated with MGR against C. ribicola, MGR in these seed families appears to be controlled by Cr4 or other R genes in very close proximity to Cr4. These associated SNPs were located in genes involved in multiple molecular mechanisms potentially underlying limber pine MGR to C. ribicola, including NBS-LRR genes for recognition of C. ribicola effectors, signaling components, and a large set of defense-responsive genes with potential functions in plant effector-triggered immunity (ETI). Interactions of associated loci were identified for MGR selection in trees with complex genetic backgrounds. SNPs with tight Cr4-linkage were further converted to TaqMan assays to confirm their effectiveness as MAS tools. This work demonstrates the successful translation and deployment of molecular genetic knowledge into specific MAS tools that can be easily applied in a selection or breeding program to efficiently screen MGR against WPBR in Alberta limber pine populations.

Characterization of Cronartium ribicola dsRNAs reveals novel members of the family Totiviridae and viral association with fungal virulence.

BACKGROUND: Mycoviruses were recently discovered in the white pine blister rust (WPBR) fungus Cronartium ribicola (J.C. Fisch.). Detection and characterization of their double stranded RNA (dsRNA) would facilitate understanding of pathogen virulence and disease pathogenesis in WPBR systems. METHODS: Full-length cDNAs were cloned from the dsRNAs purified from viral-infected C. ribicola, and their cDNA sequences were determined by DNA sequencing. Evolutionary relationships of the dsRNAs with related mycoviruses were determined by phylogenetic analysis. Dynamic distributions of the viral RNAs within samples of their fungal host C. ribicola were investigated by measurement of viral genome prevalence and viral gene expression. RESULTS: In this study we identified and characterized five novel dsRNAs from C. ribicola, designated as Cronartium ribicola totivirus 1-5 (CrTV1 to CrTV5). These dsRNA sequences encode capsid protein and RNA-dependent RNA polymerase with significant homologies to dsRNA viruses of the family Totiviridae. Phylogenetic analysis showed that the CrTVs were grouped into two distinct clades. CrTV2 through CrTV5 clustered within the genus Totivirus. CrTV1 along with a few un-assigned dsRNAs constituted a distinct phyletic clade that is genetically distant from presently known genera in the Totiviridae family, indicating that CrTV1 represents a novel genus in the Totiviridae family. The CrTVs were prevalent in fungal samples obtained from infected western white pine, whitebark pine, and limber pines. Viral RNAs were generally expressed at higher levels during in planta mycelium growth than in aeciospores and urediniospores. CrTV4 was significantly associated with C. ribicola virulent pathotype and specific C. ribicola host tree species, suggesting dsRNAs as potential tools for dissection of pathogenic mechanisms of C. ribicola and diagnosis of C. ribicola pathotypes. CONCLUSION: Phylogenetic and expression analyses of viruses in the WPBR pathogen, C. ribicola, have enchanced our understanding of virus diversity in the family Totiviridae, and provided a potential strategy to utilize pathotype-associated mycoviruses to control fungal forest diseases.

Identification and Functional Characterization of an Effector Secreted by Cronartium ribicola.

Cri-9402 was identified as a protein effector from Cronartium ribicola, based on the effect of its expression on growth of Pseudomonas syringae Psm ES4326 introduced into transiently transformed tobacco leaves and stably transformed Arabidopsis seedlings. In tobacco leaves transiently expressing its coding sequence, growth of P. syringae Psm ES4326 was inhibited. Expression of pathogenesis-related (PR) protein 2 (PR2), PR4a, endochitinase B, hypersensitive-related 201 (HSR201), HSR203J, and proteinase inhibitor 1 was upregulated but expression of PR1, coronatine insensitive 1, and abscisic acid 1 was significantly suppressed. In transformed Arabidopsis seedlings, the effector stimulated growth of P. syringae Psm ES4326; significantly suppressed expression of PR1, PR2, nonexpresser of pathogenesis-related genes 1 (NPR1), NPR3, NPR4, phytoalexin deficient 4, and salicylic acid induction deficient 2; and enhanced expression of plant defensin 1.2 (PDF1.2). The above results showed that the majority of responses to this effector in tobacco leaves were converse to those in transformed Arabidopsis. We could conclude that Cri-9402 promoted disease resistance in tobacco leaves and disease susceptibility in Arabidopsis seedlings. Its transcript was mainly expressed in aeciospores of C. ribicola and was probably involved in production or germination of aeciospores, and it was an effector potentially functioning in white pine-blister rust interactions.

Limber pine (Pinus flexilis James) genetic map constructed by exome-seq provides insight into the evolution of disease resistance and a genomic resource for genomics-based breeding.

Limber pine (Pinus flexilis) is a keystone species of high-elevation forest ecosystems of western North America, but some parts of the geographic range have high infection and mortality from the non-native white pine blister rust caused byCronartium ribicola. Genetic maps can provide essential knowledge for understanding genetic disease resistance as well as local adaptation to changing climates. Exome-seq was performed to construct high-density genetic maps in two seed families. Composite maps positioned 9612 unigenes across 12 linkage groups (LGs). Syntenic analysis of genome structure revealed that the majority of orthologs were positional orthologous genes (POGs) with localization on homologousLGs among conifer species. Gene ontology (GO) enrichment analysis showed relatively fewer constraints forPOGs with putative roles in adaptation to environments and relatively more conservation forPOGs with roles in basic cell function and maintenance. The mapped genes included 639 nucleotide-binding site leucine-rich repeat genes (NBS-LRRs), 290 receptor-like protein kinase genes (RLKs), and 1014 genes with potential roles in the defense response and induced systemic resistance to attack by pathogens. Orthologous loci for resistance to rust pathogens were identified and were co-positioned with multiple members of theR gene family, revealing the evolutionary pressure acting upon them. This high-density genetic map provides a genomic resource and practical tool for breeding and genetic conservation programs, with applications in genome-wide association studies (GWASs), the characterization of functional genes underlying complex traits, and the sequencing and assembly of the full-length genomes of limber pine and relatedPinus species.

Characterization of Heterobasidion occidentale transcriptomes reveals candidate genes and DNA polymorphisms for virulence variations.

Characterization of genes involved in differentiation of pathogen species and isolates with variations of virulence traits provides valuable information to control tree diseases for meeting the challenges of sustainable forest health and phytosanitary trade issues. Lack of genetic knowledge and genomic resources hinders novel gene discovery, molecular mechanism studies and development of diagnostic tools in the management of forest pathogens. Here, we report on transcriptome profiling of Heterobasidion occidentale isolates with contrasting virulence levels. Comparative transcriptomic analysis identified orthologous groups exclusive to H. occidentale and its isolates, revealing biological processes involved in the differentiation of isolates. Further bioinformatics analyses identified an H. occidentale secretome, CYPome and other candidate effectors, from which genes with species- and isolate-specific expression were characterized. A large proportion of differentially expressed genes were revealed to have putative activities as cell wall modification enzymes and transcription factors, suggesting their potential roles in virulence and fungal pathogenesis. Next, large numbers of simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs) were detected, including more than 14 000 interisolate non-synonymous SNPs. These polymorphic loci and species/isolate-specific genes may contribute to virulence variations and provide ideal DNA markers for development of diagnostic tools and investigation of genetic diversity.

Profiling methyl jasmonate-responsive transcriptome for understanding induced systemic resistance in whitebark pine (Pinus albicaulis).

KEY MESSAGE: RNA-seq analysis on whitebark pine needles demonstrated that methyl jasmonate (MeJA)-triggered transcriptome re-programming substantially overlapped with defense responses against insects and fungal pathogens in Pinus species, increasing current knowledge regarding induced systemic resistance (ISR) to pathogens and pests in whitebark pine. Many whitebark pine populations are in steep decline due to high susceptibility to mountain pine beetle and the non-native white pine blister rust (WPBR). Resistance, including induced systemic resistance (ISR), is not well characterized in whitebark pine, narrowing the current options for increasing the success of restoration and breeding programs. Exogenous jasmonates are known to trigger ISR by activating the plant's immune system through regulation of gene expression to produce chemical defense compounds. This study reports profiles of whitebark pine needle transcriptomes, following methyl jasmonate (MeJA) treatment using RNA-seq. A MeJA-responsive transcriptome was de novo assembled and transcriptome profiling identified a set of differentially expressed genes (DEGs), revealing 1422 up- and 999 down-regulated transcripts with at least twofold change (FDR corrected p < 0.05) in needle tissues in response to MeJA application. GO analysis revealed that these DEGs have putative functions in plant defense signalling, transcription regulation, biosyntheses of secondary metabolites, and other biological processes. Lineage-specific expression of defense-related genes was characterized through comparison with MeJA signalling in model plants. In particular, MeJA-triggered transcriptome re-programming substantially overlapped with defense responses against WPBR and insects in related Pinus species, suggesting that MeJA may be used to improve whitebark pine resistance to pathogens/pests. Our study provides new insights into molecular mechanisms and metabolic pathways involved in whitebark pine ISR. DEGs identified in this study can be used as candidates to facilitate identification of genomic variation contributing to host resistance and aid in breeding selection of elite genotypes with better adaptive fitness to environmental stressors in this endangered tree species.

Saturated genic SNP mapping identified functional candidates and selection tools for the Pinus monticola Cr2 locus controlling resistance to white pine blister rust.

Molecular breeding incorporates efficient tools to increase rust resistance in five-needle pines. Susceptibility of native five-needle pines to white pine blister rust (WPBR), caused by the non-native invasive fungus Cronartium ribicola (J.C. Fisch.), has significantly reduced wild populations of these conifers in North America. Major resistance (R) genes against specific avirulent pathotypes have been found in several five-needle pine species. In this study, we screened genic SNP markers by comparative transcriptome and genetic association analyses and constructed saturated linkage maps for the western white pine (Pinus monticola) R locus (Cr2). Phenotypic segregation was measured by a hypersensitive reaction (HR)-like response on the needles and disease symptoms of cankered stems post inoculation by the C. ribicola avcr2 race. SNP genotypes were determined by HRM- and TaqMan-based SNP genotyping. Saturated maps of the Cr2-linkage group (LG) were constructed in three seed families using a total of 34 SNP markers within 21 unique genes. Cr2 was consistently flanked by contig_2142 (encoding a ruvb-like protein) and contig_3772 (encoding a delta-fatty acid desaturase) across the three seed families. Cr2 was anchored to the Pinus consensus LG-1, which differs from LGs where other R loci of Pinus species were mapped. GO annotation identified a set of NBS-LRR and other resistance-related genes as R candidates in the Cr2 region. Association of one nonsynonymous SNP locus of an NBS-LRR gene with Cr2-mediated phenotypes provides a valuable tool for marker-assisted selection (MAS), which will shorten the breeding cycle of resistance screening and aid in the restoration of WPBR-disturbed forest ecosystems.

Genetic Diversity and Population Structure of Whitebark Pine (Pinus albicaulis Engelm.) in Western North America.

Whitebark pine (WBP, Pinus albicaulis Engelm.) is an endangered conifer species due to heavy mortality from white pine blister rust (WPBR, caused by Cronartium ribicola) and mountain pine beetle (Dendroctonus ponderosae). Information about genetic diversity and population structure is of fundamental importance for its conservation and restoration. However, current knowledge on the genetic constitution and genomic variation is still limited for WBP. In this study, an integrated genomics approach was applied to characterize seed collections from WBP breeding programs in western North America. RNA-seq analysis was used for de novo assembly of the WBP needle transcriptome, which contains 97,447 protein-coding transcripts. Within the transcriptome, single nucleotide polymorphisms (SNPs) were discovered, and more than 22,000 of them were non-synonymous SNPs (ns-SNPs). Following the annotation of genes with ns-SNPs, 216 ns-SNPs within candidate genes with putative functions in disease resistance and plant defense were selected to design SNP arrays for high-throughput genotyping. Among these SNP loci, 71 were highly polymorphic, with sufficient variation to identify a unique genotype for each of the 371 individuals originating from British Columbia (Canada), Oregon and Washington (USA). A clear genetic differentiation was evident among seed families. Analyses of genetic spatial patterns revealed varying degrees of diversity and the existence of several genetic subgroups in the WBP breeding populations. Genetic components were associated with geographic variables and phenotypic rating of WPBR disease severity across landscapes, which may facilitate further identification of WBP genotypes and gene alleles contributing to local adaptation and quantitative resistance to WPBR. The WBP genomic resources developed here provide an invaluable tool for further studies and for exploitation and utilization of the genetic diversity preserved within this endangered conifer and other five-needle pines.

Genetic mapping of Pinus flexilis major gene (Cr4) for resistance to white pine blister rust using transcriptome-based SNP genotyping.

BACKGROUND: Linkage of DNA markers with phenotypic traits provides essential information to dissect clustered genes with potential phenotypic contributions in a target genome region. Pinus flexilis E. James (limber pine) is a keystone five-needle pine species in mountain-top ecosystems of North America. White pine blister rust (WPBR), caused by a non-native fungal pathogen Cronartium ribicola (J.C. Fisch.), has resulted in mortality in this conifer species and is still spreading through the distribution. The objective of this research was to develop P. flexilis transcriptome-wide single nucleotide polymorphism (SNP) markers using RNA-seq analysis for genetic mapping of the major gene (Cr4) that confers complete resistance to C. ribicola. RESULTS: Needle tissues of one resistant and two susceptible seedling families were subjected to RNA-seq analysis. In silico SNP markers were uncovered by mapping the RNA-seq reads back to the de novo assembled transcriptomes. A total of 110,573 in silico SNPs and 2,870 indels were identified with an average of 3.7 SNPs per Kb. These SNPs were distributed in 17,041 unigenes. Of these polymorphic P. flexilis unigenes, 6,584 were highly conserved as compared to the genome sequence of P. taeda L (loblolly pine). High-throughput genotyping arrays were designed and were used to search for Cr4-linked genic SNPs in megagametophyte populations of four maternal trees by haploid-segregation analysis. A total of 32 SNP markers in 25 genes were localized on the Cr4 linkage group (LG). Syntenic relationships of this Cr4-LG map with the model conifer species P. taeda anchored Cr4 on Pinus consensus LG8, indicating that R genes against C. ribicola have evolved independently in different five-needle pines. Functional genes close to Cr4 were annotated and their potential roles in Cr4-mediated resistance were further discussed. CONCLUSIONS: We demonstrated a very effective, low-cost strategy for developing a SNP genetic map of a phenotypic trait of interest. SNP discovery through transcriptome comparison was integrated with high-throughput genotyping of a small set of in silico SNPs. This strategy may be applied to mapping any trait in non-model plant species that have complex genomes. Whole transcriptome sequencing provides a powerful tool for SNP discovery in conifers and other species with complex genomes, for which sequencing and annotation of complex genomes is still challenging. The genic SNP map for the consensus Cr4-LG may help future molecular breeding efforts by enabling both Cr4 positional characterization and selection of this gene against WPBR.

Transcriptome analysis of the white pine blister rust pathogen Cronartium ribicola: de novo assembly, expression profiling, and identification of candidate effectors.

BACKGROUND: The fungus Cronartium ribicola (Cri) is an economically and ecologically important forest pathogen that causes white pine blister rust (WPBR) disease on five-needle pines. To cause stem cankers and kill white pine trees the fungus elaborates a life cycle with five stages of spore development on five-needle pines and the alternate host Ribes plants. To increase our understanding of molecular WP-BR interactions, here we report genome-wide transcriptional profile analysis of C. ribicola using RNA-seq. RESULTS: cDNA libraries were constructed from aeciospore, urediniospore, and western white pine (Pinus monticola) tissues post Cri infection. Over 200 million RNA-seq 100-bp paired-end (PE) reads from rust fungal spores were de novo assembled and a reference transcriptome was generated with 17,880 transcripts that were expressed from 13,629 unigenes. A total of 734 unique proteins were predicted as a part of the Cri secretome from complete open reading frames (ORFs), and 41 % of them were Cronartium-specific. This study further identified a repertoire of candidate effectors and other pathogenicity determinants. Differentially expressed genes (DEGs) were identified to gain an understanding of molecular events important during the WPBR fungus life cycle by comparing Cri transcriptomes at different infection stages. Large-scale changes of in planta gene expression profiles were observed, revealing that multiple fungal biosynthetic pathways were enhanced during mycelium growth inside infected pine stem tissues. Conversely, many fungal genes that were up-regulated at the urediniospore stage appeared to be signalling components and transporters. The secreted fungal protein genes that were up-regulated in pine needle tissues during early infection were primarily associated with cell wall modifications, possibly to mask the rust pathogen from plant defenses. CONCLUSION: This comprehensive transcriptome profiling substantially improves our current understanding of molecular WP-BR interactions. The repertoire of candidate effectors and other putative pathogenicity determinants identified here are valuable for future functional analysis of Cri virulence and pathogenicity.

Anti-microbial peptide (AMP): nucleotide variation, gene expression, and host resistance in the white pine blister rust (WPBR) pathosystem.

Pinus monticola antimicrobial peptide (PmAMP1) inhibits growth of Cronartium ribicola and other fungal pathogens. C. ribicola causes white pine blister rust and has resulted in a dramatic reduction of native white pines across North America. Quantitative disease resistance (QDR) is a highly desirable trait screened in breeding programs for durable resistance against C. ribicola. Along with phenotyping on a collection of germplasms, we analyzed PmAMP1 transcript and protein expression and re-sequenced the full-length gene including its promoter region. A mixed linear model was used to identify the association of single nucleotide polymorphisms (SNPs) with accumulated protein and stem QDR levels. Among 16 PmAMP1 SNPs identified in the present study, we found an association of protein levels with 6 SNPs (P < 0.05), including 2 in the 5'-untranslated region (UTR), 3 in the open reading frame (ORF) region with 2 nonsynonymous SNPs, and 1 SNP in the 3'-UTR. Another set of six SNPs was associated with stem QDR levels (P < 0.05), with one localized in the promoter region and the other five in the ORF region with four nonsynonymous changes, suggesting that multiple isoforms may have antifungal activity to differing degrees. Of three common PmAMP1 haplotypes, the trees with haplotype 2 showed high QDR levels with moderate protein abundance while those trees with haplotype 3 exhibited low QDR levels in the susceptible range and the lowest level of protein accumulation. Thus, an association of gene variations with protein abundance and resistance-related traits may facilitate elucidation of physiological contribution of PmAMP1 to host resistance.

Comparative proteomic profiles of Pinus monticola needles during early compatible and incompatible interactions with Cronartium ribicola.

The proteomic profiles of primary needles from Cr2-resistant and cr2-susceptible Pinus monticola seedlings were analysed post Cronartium ribicola inoculation by 2-DE. One hundred-and-five protein spots exhibiting significant differential expression were identified using LC-MS/MS. Functional classification showed that the most numerous proteins are involved in defence signalling, oxidative burst, metabolic pathways, and other physiological processes. Our results revealed that differential expression of proteins in response to C. ribicola inoculation was genotype- and infection-stage dependent. Responsive proteins in resistant seedlings with incompatible white pine blister rust (WPBR) interaction included such well-characterized proteins as heat shock proteins (HSPs), reactive oxygen species (ROS) scavenging enzymes, and intermediate factors functioning in the signal transduction pathways triggered by well-known plant R genes, as well as new candidates in plant defence like sugar epimerase, GTP-binding proteins, and chloroplastic ribonucleoproteins. Fewer proteins were regulated in susceptible seedlings; most of them were in common with resistant seedlings and related to photosynthesis among others. Quantitative RT-PCR analysis confirmed HSP- and ROS-related genes played an important role in host defence in response to C. ribicola infection. To the best of our knowledge, this is the first comparative proteomics study on WPBR interactions at the early stages of host defence, which provides a reference proteomic profile for other five-needle pines as well as resistance candidates for further understanding of host resistance in the WPBR pathosystem.

Antifungal activity of a Pinus monticola antimicrobial peptide 1 (Pm-AMP1) and its accumulation in western white pine infected with Cronartium ribicola.

Pinus monticola antimicrobial peptide 1 (Pm-AMP1) was expressed and purified from bacterial cell lysate and its identity and purity confirmed by Western blot analysis using the Pm-AMP1 antibody. Application of Pm-AMP1 resulted in visible hyphal growth inhibition of Cronartium ribicola , Phellinus sulphurascens , Ophiostoma montium , and Ophiostoma clavigerum 3-12 days post-treatment. Pm-AMP1 also inhibited spore germination of several other phytopathogenic fungi by 32%-84% 5 days post-treatment. Microscopic examination of C. ribicola hyphae in contact with Pm-AMP1 showed distinct morphological changes. Seven western white pine ( Pinus monticola Douglas ex D. Don) families (Nos. 1, 2, 5, 6, 7, 8, 10) showing partial resistance to C. ribicola in the form of bark reaction (BR) were assessed by Western immunoblot for associations between Pm-AMP1 accumulation and family, phenotype, canker number, and virulence of C. ribicola. There was a significant difference (p < 0.001) in mean Pm-AMP1 protein accumulation between families, with higher levels detected in the full-sib BR families (Nos. 1, 2, 5) than the half-sib BR families (Nos. 6, 7). Family 8, previously described as a Mechanism 'X' BR family, had the highest number of BR seedlings and displayed high Pm-AMP1 levels, whereas the susceptible family (No. 10) showed the lowest levels (p < 0.05). Family 1 showed a significant association between Pm-AMP1 accumulation and overall seedling health (p < 0.01, R = 0.533), with higher protein levels observed in healthy versus severely infected seedlings. In general, low Pm-AMP1 levels were observed with an increase in the number of cankers per seedling (p < 0.05), and seedlings inoculated with the avirulent source of C. ribicola showed significantly higher Pm-AMP1 levels (p < 0.05) in the majority of BR families. Cis-acting regulatory elements, such as CCAAT binding factors, and an AG-motif binding protein were identified in the Pm-AMP1 promoter region. Multiple polymorphic sites were identified within the 5' untranslated region and promoter regions. Our results suggest that Pm-AMP1 is involved in the western white pine defense response to fungal infection, as observed by its antifungal activity on C. ribicola and a range of phytopathogens as well as through its association with different indicators of resistance to C. ribicola.